arXiv:1502.02883vl [cond-mat.stat-mech] 10 Feb 2015 


The Ising chain constrained to an even or odd number of positive 
spins 


Michael T. Gastner 1,2 

1 Institute of Technical Physics and Materials Science, Research Centre for Natural Sciences, Hungarian 
Academy of Sciences, P.O. Box 49, H-1525 Budapest, Hungary 

2 Department of Engineering Mathematics, University of Bristol, Merchant Venturers Building, Woodland 
Road, Bristol BS8 1UB, United Kingdom. 

E-mail: m.gastner@bristol.ac.uk 


Abstract. We investigate the statistical mechanics of the periodic one-dimensional Ising chain when the 
number of positive spins is constrained to be either an even or an odd number. We calculate the partition 
function using a generalization of the transfer matrix method. On this basis, we derive the exact magneti¬ 
zation, susceptibility, internal energy, heat capacity and correlation function. We show that in general the 
constraints substantially slow down convergence to the thermodynamic limit. By taking the thermodynamic 
limit together with the limit of zero temperature and zero magnetic field, the constraints lead to new scaling 
functions and different probability distributions for the magnetization. We demonstrate how these results 
solve a stochastic version of the one-dimensional voter model. 


1. Introduction 

For almost one century, the Ising model of ferromagnetism has been a cornerstone of statistical mechanics [T]. 
It is one of very few problems that can, at least in one and two dimensions, be solved exactly [2|. Its 
applications range from solid state physics [3] over neuroscience j4] to collective social phenomena [5] . In its 
basic form the Ising model is based on the Hamiltonian 

N N 

E(cr) = -J^(7 i cr i+ i - H^ai, (1) 

i= 1 z=l 

where each spin in the vector er = (ui,..., cr/v) can take only the values ±1 and we assume periodic boundary 
conditions so that cjn+i = txi- The parameter J is the strength of interactions between spins and H an 
external magnetic field. We allow J to take positive and negative values, thereby considering both the ferro- 
and antiferromagnetic case. 

Many generalizations of the model have been investigated since Ising’s groundbreaking publication, for 
example long-range interactions [5j, spin glasses [?i and permitting more than two possible spin states [S]. 
In this article we investigate two different variations of the Ising model. First, we restrict the number of 
positive spins 

V + (<r) = i|v + X>^ (2) 

to an even number. That is, the Hamiltonian is given by Eq. [I] if N + is even and E = oo if it is odd. In the 
second model, Eq. [T] holds if N + is odd, whereas an even N + is forbidden. We will refer to these two models 
as “even” or “odd” Ising model respectively 0. 

J Sometimes the term “odd Ising model” is used for a spin glass model by Villain m which is unrelated to our work. 


In Sec. mm we will motivate the even and odd model by showing that they are equivalent to a simple 
opinion formation model. In Sec. [5] we demonstrate how the transfer matrix method for the unconstrained 
Ising model can be modified to derive the partition functions of the even and odd model. Section [6] contains 
a derivation of the magnetization and susceptibility of both models. We deduce the nearest-neighbour 
correlations, internal energy and heat capacity in Sec. [7] and the correlation function in Sec. [S] As we show 
in Sec. I9land [l0l the constrained models approach the thermodynamic limit in a different manner than the 
usual unconstrained model when the temperature and magnetic field simultaneously go to zero. We apply 
these results to the opinion formation model in Sec. [TT] before summarizing the key findings in Sec. 1121 

Before proceeding, we emphasize that N + is not a fixed number, neither in the even nor odd model. It is 
still permitted to take a multitude of values (e.g. in the even model N + = 0, 2,4,..., 2|_iV/2j), but with the 
restriction that configurations with either odd or even N + are excluded. In a Monte Carlo simulation, this 
restriction could be imposed by initializing the spins with an even or odd N + and subsequently flipping two 
distinct spins simultaneously in each update. Because such a Markov chain is not ergodically exploring the 
configurations of the conventional (i.e. unconstrained) Ising model, we should not expect that the equilibrium 
properties are equal. One purpose of this article is to convince ourselves that the thermodynamic limits (i.e. 
N —> oo) of the even and odd models are indeed the limits of the unconstrained model for fixed temperature. 
However, we will point out differences when the thermodynamic limit is taken simultaneously with the limit 
of zero temperature and zero magnetic field. 

2. Motivation: Stochastic synchronous voter model 

We consider a version of the voter model with stochastic opinion updates. Individuals are placed on the 
N sites of a one-dimensional chain with periodic boundary conditions. Each individual holds one of two 
possible opinions: “black” or “white”. We associate each site i with a binary variable Wi whose values are 

I 1 if z is black at time t , 

Wi(t) = (3) 

1—1 if z is white at time t. 

At each discrete time step t, all individuals synchronously update their opinions 0. (We will discuss 
asynchronous updates in Sec. [4]) Each individual randomly chooses one of their two nearest neighbours and 
adopts her opinion with probability or chooses the opposite opinion with probability p- = 1 — p+. Thus, 
the probability that i’s next opinion is Q = ±1 can be expressed as 

f P+ if § (<*>i-i(t) +Ui+i(t)) = 1, 

PrMi + 1) = fi|wi_i(t),w i+ i(i)] = l \ if Wi_i(t) + u l+1 (t) = 0, (4) 

[p~ if § +w i+ i(t)) = -1, 

where the subindices are interpreted modulo N to satisfy the periodic boundary conditions. 

What are the equilibrium properties of this model? For example, how many pairs of neighbours will on 
average disagree? And what are the typical fluctuations around this average value? We will demonstrate 
that these questions can be analytically answered by mapping the problem to an Ising model on the dual 
lattice with an even number of negative spins. (We will explain the origin of the even-numbered constraint 
in Sec. [3]) For even (odd) N, the opinion model will consequently map onto the even (odd) Ising model. 

Let us first clarify that the variables ui t cannot directly be interpreted as Ising spins a t . For simplicity’s 
sake, let us assume for a moment that N is even. In the limiting case of p+ = 1, there are two stationary 
states where all sites have reached either a black or white consensus (Fig. QJi and [[):>). For synchronous 
updates there is, however, also a periodic state where the opinions alternate in space m-- if all odd sites are 
black and all even sites white at time i, all opinions are inverted at t + 1 and return to the original state at 
t + 2 (Fig. [TJ;). Unlike in the zero-temperature Ising model, we thus have apparently more than two ground 
states. 
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Figure 1. The (a), (b) stationary and (c) periodic states of the opinion dynamics of Eq. [T]in the limit 
p .|_ = 1 with synchronous updates. A site i with iv; = ±1 is represented by a black (white) square. The 
spins of the associated Ising model are shown as + or — signs above the links. 


We can, however, establish a connection to the Ising model if we assign spins cq to the i-th link (i.e. 
between the sites i and i + 1) rather than the sites themselves. We set tr* = 1 if both sites connected by the 
link agree and a t = —1 if they disagree, 

<Ji = uwj+i. (5) 

In terms of cr,, both consensus states are mapped to maximally positive magnetization, whereas in the 
alternating state all spins are negative. Thus, the limit p+ = 1 can be mapped to the zero-temperature 
ferromagnetic Ising model. We will now argue that for any 0 < < 1, there is a finite-temperature Ising 

model whose equilibrium properties are those of the original opinion dynamics given by Eq. [I] 


3. Mapping the synchronous voter model to an odd or even Ising model 


wjy . What is the probability Pr(A —► B ) to find the 


Suppose that the opinions at time t are uq 
opinions ,..., at time t + 1? Assuming that the probabilities in Eq. [J] are independent for all *, 


N 


Pr(A B) = Pr 


,(B) 


,.(A) (A) 

1) U i +1 


We want to show that 
Pr(A - 


B) _ p(e ( - a ’>-e ( - b '>) 


Pr (B —> A) 

where E^ is the energy of the spins cr- A ^ = in the Ising model without magnetic field, 


£(^) =-jfv (A V (A) 


i =1 


i+ 1> 


and similarly for state B. Furthermore, 
In (2 y/p^pl) 


P = — 


2 J 


( 6 ) 


( 7 ) 


( 8 ) 


(9) 


so that every p+ can be mapped to a temperature ( ksP ) , where is the Boltzmann constant. Equation!?] 

is the detailed balance condition for the Ising model m- Consequently, the equilibrium properties of the 
spins cti can be deduced from the model’s partition function. 

Before deriving Eq. [T1 we emphasize that not all spin configurations are possible. The number of negative 
spins must be even; otherwise the opinions Wj would change an odd number of times as we go once through 
the chain so that we would not end up with the same opinion with which we started. The restriction to an 


3 



















even number of negative spins changes the partition function of this model compared to the unconstrained 
Ising model. 

First, however, we still need to justify Eq. [7] Let us denote the number of neighbouring spins with 
opposite signs in states A and B by n^ and n^ respectively. Because E^ = J (2n < -' 4 ) — iV) and 
E = J (2— ./V), we can rewrite the right-hand side of Eq. [7] as 


e /3(_E (yt) -£; (fl) ) _ p 2/3J(ra (A) -rt (B) ) 

( 10 ) 

Because of Eq. EJandEl only factors p + , \ and p- can appear in Pr(A —► B) and Pr(H —>• A), 

Pr (A^ B) =p a + 1 p a _ 2 /2 a3 , 

(ii) 

Pr (B -> A) =p b {p b */2 b3 . 

( 12 ) 

Assuming p + ^ 1, the exponents cq, bi are uniquely determined, 
correct.) There is one factor for each site, so 

(If p+ = 1 and thus (3 = 0, Eq.[7]is trivially 

a\ + a 2 + 03 = 61 + b 2 + b 3 = N. 

(13) 

Because OjCj+i = —1 if and only if w* + u>i+ 2 = 0, Eq. [4] implies 


a 3 = n (A) , 

(14) 

63 = rS B \ 

(15) 

From Eq. [4] it also follows that 


„ „ _ 1 VW ,(B) ( (A) (A)\ 

a\ - a 2 — - 2_, (^A-i + w *+ 1 J i 

(16) 

h h - 1 - , (j4) (, 

(17) 

Splitting the sums and shifting the summation index shows that the sums in Eq. ITOl and H71 are equal, thus 

a\ — a 2 = bi — b 2 . 

(18) 

Combining Eq. [T31 El E] and El 


a 2 = N — a\ — rS A \ 

(19) 

bi = a\ + i (n ( ' 4) - , 

( 20 ) 

b 2 = N — gh — i ^ n ^ + n^ B ^ , 

( 21 ) 

so that, by plugging into Eq. [Til and fl2l we obtain 


Pr(A —y B) , _ \n^ B ^—n^ A ^ 

Pl(B A) - 

(22) 


Comparing Eq. El [TO] and [22] proves Eq. [7] 


4. The voter model with random asynchronous updates 

Not only the voter model with perfectly synchronous updates of opinions can be mapped to an even or odd 
Ising model. We will now argue that, by defining the spins as in Eq. [5j we can also interpret asynchronous 
updates of randomly selected single opinions in terms of an Ising Hamiltonian. While the synchronous case, 
as shown in the previous section, corresponds to a positive spin interaction J and zero magnetic field H , the 
asynchronous case leads to J = 0 and, in general, H ^ 0 for the following reason. 

Suppose opinion cu, is chosen to be updated. The probability to have opinion Q in the next time step is 
given by Eq. [4] while all other opinions remain unchanged. Then the only spins affected are cq_ 1 and cq so 
that we can ignore the rest of the chain. If uii changes between states A and I?, then we can distinguish the 
three cases depicted in Fig. [2] either 


4 




Pr = 1/2 



Pr = 1/2 


(ii) 


(iii) 


P+ 



P+ 


Figure 2. Transition probabilities for asynchronous updates. Depicted are three representative cases 
where only opinion i changes between states A and B. All other cases can be generated by inverting all 
opinions (from white to black and vice versa) and/or interchanging the order of the chain so that i — 1 and 
* + 1 trade places. 


(i) Pr (A -a B) = Pr(B ->A) = i and a[ A) + a (A \ = a {B) + a™ = 0 or 

(ii) Pr(A —» B) = p + and Pr (B -A A) = p_ and a— a^ B \ = —2 or 

(iii) Pr(A —/ B) = p- and Pr(B — > A) = p + and + a^ A \ = = 2. 

In summary, we can write all of these cases as 


Pr(A -A B) 
Pr (B -A- A) 




which is of the form of Eq. [7] with the energy E = —H er; and inverse temperature 


(23) 


kAAeA Ez>. (24) 

If p + £ (i, 1), we must have H > 0 to obtain a positive temperature. Asynchronous opinion updates 
then tend to favour the states depicted in Fig. Oja) an d (b) where the spins are all positive. Generally, the 
states in Fig. He) are suppressed whenp + >1/2 and the dynamic rule mixes synchronous and asynchronous 
updates (e.g. by updating a fraction of the opinions in every update as in Ref. mi)- The opposite is true for 
p. <1/2 where asynchronous updates generate sequences of alternating opinions and suppress unanimity. 
All cases, however, have in common that the periodic boundary conditions in the opinions generate an even 
number of negative spins, resulting in an even (odd) Ising model for even (odd) N. 

Whether synchronous, asynchronous or partially synchronous updates are more realistic depends on the 
situation one wishes to model. Asynchronous updates have a long tradition in physics (e.g. the Glauber or 
Metropolis rules for dynamic Ising models), but synchronous updates, especially in the context of stochastic 
cellular automata HE have also been investigated (for example in [141 IT51ITB1 ). If agents can only make 
decisions at discrete times (e.g. only at the end of a business day or if biological populations exhibit strongly 
peaked cyclic activity mi), then synchronous or partially synchronous updates are more applicable. Here 
we do not intend to argue for any particular update rule. Generally, one has to be humble about social or 
economic interpretations of such simple rules m because true opinion dynamics is far more complex. Our 
focus here is rather on the model’s structural properties in order to motivate how the even and odd Ising 
models can arise from another two-state model. 
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5. Partition function 


We denote the partition function for the even and odd Ising model by Z e and Z 0 respectively, 


Z e = Y 

cr with 
even 7 V_|.(<t) 

Z 0 = Y e- 0EM . 

cr with 
odd iV_|_(<r) 


(25) 

(26) 


If we associate a spin <jj = 1 with the bra vector (+1| = (1,0) and <r» = —1 with (—1| = (0,1), we can write 
Z e with the transfer matrix of the unconstrained model nn 


/ e {t{J+H) e -/3 J \ 

^ e -PJ e P{J-H) J 


as 




Y (oil P |o 2 ) <cr 2 | p 1 0 - 3 ) • ■ ■ (owl p |cti) = Tr 

cr with 
even -/V_|_(<x) 


( \ 

Y |cti)(oi|P...|ctjv)(o-w|P 

cr with 

\even iV_|_(<x) J 


(27) 


(28) 


Similarly, 


Let us define 


Z Q = Tr 


( \ 

Y |oi) (oi| P ... |ow) (ow| P 

cr with 

\odd AT_|_(<r) / 


Me = Y |<7l) (<7l| P . . . |tTjv) (CTJV| P, 

cr with 
even iV_|_(cr) 

M 0 = Yj |oi) (cti| P ... |tJw) (ow| P- 

cr with 
odd iV_|_(<x) 


Induction on N proves 


( M e , M 0 \ / 1-1) (-1| P, | + 1)(+1|P\ 

V Mo, m £ j Vi+i)(+i|p. 1-1) (-11 py 


With the definition 


( hi) (-11 P, l+l)(+l|P 
V I+1)(+1|P. 1-1) (-11P 


0 

0 

e 0 (J+H) 

e- 8J \ 

e -/ 3 J 

e P (J-H) 

0 

0 

P (J+H) 

e -P J 

0 

0 

0 

0 

e -P J 

e p (j-h) y 


we can write 

Ze = Tr (M e ) = i Tr (Q^) ■ 
To simplify the notation further, we introduce 


(29) 

(30) 

(31) 

(32) 

(33) 

(34) 

(35) 


x = pH, 

y=pj. 
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(36) 

(37) 








The eigenvalues of Q are then 

Ai ,2 = e y ^cosh a; ± \J sinh 2 x + e~ 4y ^j , (38) 

•^ 3,4 = e v sinh a; ± \J cosh 2 x — e~ iy ^j . (39) 

Consequently, 

Z e = - (Af + Af + Af + Af) . (40) 

We can derive Z 0 as follows. The eigenvalues Aj .,2 are also the eigenvalues of P and therefore the partition 
function of the unconstrained Ising model is = Af + Af. Moreover Z u = Z e + Z a so that 

= \ (Af + Af - Af - Af) . (41) 

Because Ai is the leading eigenvalue, we find in the thermodynamic limit (i.e. N —> oo) with fixed x and y 
that Z e oc Z 0 oc Z u oc Af. As a consequence, all equilibrium properties of the even and odd Ising models 
converge to the same limits as the unconstrained model. However, we will analytically derive in Sec. [9] 
different scaling limits for N —> oo when temperature and magnetic field go to their critical value (i.e. zero) 
such that Ne~ 2y and N sinh x are asymptotically constants. For this purpose, it will be instructive to derive 
first some exact formulae for finite N. 


6. Magnetization and Susceptibility 


We first calculate the mean magnetization per spin 

(42) 

where Z is the partition function of the model in question and the angle brackets denote the ensemble 
average. With the auxiliary functions 

sinh x 


si(x,y) = 


ci(x,y) = 

we can write Eq. 25] as 


( m) ei0 (x,y ) = 


\J sinh 2 x + e~ 4y 
cosher 

\/ cosh 2 x — e~ Ay 

si (Af-A^Td (Af-Af) 


(43) 

(44) 

(45) 


Af + Af ± A^ ± Af 

where the upper signs apply to the even and the lower signs to the odd model. 

In the special case x = 0, applicable to the synchronous voter model, we insert the eigenvalues from 
Eq. [55] and [55] (Fig. [5J.i) 


(m} ei0 {x = 0 , 2 /) = 


(' m) u {x = 0,y) = 0 


T; 


' sinh(2y) ^ 


(JV-l)/2 


if N is even, 
if N is odd. 


(46) 


cosh A 7/+sinh /v y \ 2 J 

Hence, for odd N, even when there is no external magnetic field (i.e. H = 0), the magnetization is generally 
different from zero. This phenomenon can be intuitively explained. The constraint of an even number of 
positive spins prevents for odd N a ground state with perfectly aligned positive spins. However, the state 
with <j i =... = ctjv = —1 is permitted and therefore the mean magnetization in the limit y —> oo is — 1. The 
same argument applies with opposite signs to the odd model. In the antiferromagnetic limit (i.e. y —> —oo) 
the neighbouring spins prefer to be in opposite directions, but an odd N forces at least one pair to point in 
the same direction and thus m ei0 = (—1 )^ iv±1 )/ 2 /7V. 
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Figure 3. The mean magnetization (ra e ) (a) as a function of y when x = 0, (b) as a function of x when 
y = 0. (c) The mean magnetization (m 0 ) for y — 0. 


The relevant case for the asynchronous voter model is y = 0 where the interactions between spins are 
negligible compared to the external magnetic field, 


(m) e ,o(x 1 y = 0) = 


sinh(2:r) coslr v 2 x ± (—l) w sinh i ' 


N 


(47) 


2 cosn x ± (— 1)^ sinh 1 ’' x 

The functions are plotted in Fig. [3])b) and (c). In the unconstrained model the magnetization ( m) u (x,y = 
0) = tanhs is independent of N. However, in the even and odd models, the constraints on the number 
of spins acts as an effective interaction so that the partition function does not factorize although J = 0. 
Consequently, the magnetization of Eq. [47] depends on N. 

The fluctuations in the magnetization are measured by the susceptibility per spin 

d ( to ) 


X = /3N ((m 2 ) - (m) 2 ) = 


dH 


(48) 


Taking the derivative of Eq. [?S] for general x and y is in principle possible, but leads to rather lengthy 
expressions. We focus here instead directly on the two special cases x = 0 and y = 0. 

For x = 0 (i.e. in the absence of an external magnetic field), 

' /3e 2y (cosh JV y-sinh" y±2~ N/2 N sinh w / 2 “ 1 (2y)) 

(coslv v / 2 y±sinh ;V / 2 yj 2 

/3e 2y (cosh 2iv y- sinh 2iv y-2 1 ~ N Nsinh N ~ 1 (2y)) 


Xe,o( x = 0 ,y) = 


(cosh^ y+sinh^ y) 


if N is even, 
if N is odd, 


(49) 


compared to Xu = /3e 2y ( cosh w y — sinh ,v y)/ (cosh jV y + sinh" v y). Plotting Xe,o in Fig. [4|i and [4}o, the most 
striking feature for odd N is limy^oo Xe,o = 0, whereas the unconstrained Ising model (and the even model 
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(a) 6 



Figure 4. The susceptibility x divided by the inverse temperature /3 for (a), (b) zero magnetic field H as 
a function of y = fi.J and (c), (d) zero spin interaction J as a function of x = 0H. Panels (a) and (c) show 
the results for the even Ising model, (b) and (d) for the odd model. 


for even N) reaches in this limit its maximum susceptibility (3N. The reason is that, as already mentioned, 
the even and odd models for odd N only have one ground state each, but the unconstrained model has two. 

For the odd model with even N we observe yet another interesting phenomenon. The susceptibility 
reaches its maximum in the limit y —> oo, but with a smaller value than the unconstrained or even models, 
namely liniy-xx, \o = /3(N 2 — 4)/(3 N). The explanation is that the perfectly aligned ground states of the 
unconstrained models are not permitted. Therefore, the states of minimum energy in the odd model are the 
first excited states of the unconstrained model whose magnetization is not confined to the extreme values 
m = ± 1 . 

In the case of no internal interactions (i.e. y = 0), 

/ m 0 /^osh^zTt-l^sinh^z . {-l) N Nsmh N ~ 2 {2x) 

Xe,o(x, y = 0) = /? - - ± --- -2 

^ cosh x ± (—1)-™ smh x 2^- 2 (cosh JV x±(-l)^sinh Ar x) 

while Xu = /?cosh -2 x. We plot \e,o in Fig.HJ: and0J:l. For odd N, they satisfy Xe(x, 0) = X 0 (—x, 0) because 
in this case the odd model is equivalent to the even model with flipped signs of spins and magnetic field. 
If N is even, Xe an d Xo ar e intrinsically symmetric, but with larger values in the tails of the even model 
because Xe/Xo 3 N/(N — 2) as |x| —► oo and y = 0. 
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7. Nearest-neighbour correlations, internal energy and heat capacity 


If we replace in Eq. [42]the partial derivative with respect to x by differentiation with respect to y , we obtain 
the mean nearest-neighbour correlation 



If we define the functions 

s 2 {x 1 y) = 2e~ 3y (sinh 2 x + e ~ 4y ) , (52) 

c 2 (x,y) = 2e~ 3y (cosh 2 x — e~ 4y ) ^ , (53) 


then 


(gi)e,o( x >y) = 1 


„ (\N—1 \N -IS I t\N-l \N-1\ 

s 2 (A 2 - X 1 ) ± c 2 (A 3 - A 4 ) 


Without external magnetic field, 

(gi)e,o( x = y) = 


Af + A^ ± Af ± Af 


-i _ cosh^ 2 1 y^sinh-^/ 2 1 y 
ey (cosh^/ 2 yisinh^/ 2 y'j 

Wu{x = 0,y) = l + ^v^ 


if N is even, 
if N is odd. 


(54) 


(55) 


If N is odd, ( gi) et0 (x = 0 ,y) is equal to the correlation in the unconstrained model for the following reason. 
The spin configurations in the unconstrained model can be divided into two sets: one set containing all 
configurations of the even model and another set with all odd-numbered states. We can map every element 
in one set uniquely to the configuration in the other set that has all spins inverted. Because the sum of 
Eq. m is invariant if all spins are simultaneously flipped, the average correlations must be equal in both 
sets. The same argument cannot be applied to even N, however, because inverting the spins in the even or 
odd set generates another spin in the same set. As a consequence, (g\) e and (g i) Q are in this case different 
functions. 

With a magnetic field, but with vanishing spin interactions, 


(gi)e,o(x,y = 0 ) = 1 - 


cosh N 2 x =F (—1)^ sinh w 2 x 
cosh w x ± (—1)^ sinh^ x 


(56) 


compared to (gi) u = tanh 2 x. For odd N, we find (gi) e (x, 0) = (gi) 0 (— x, 0) for the same reason as 
discussed after the corresponding Eq. [50] for the susceptibility. We also find again that, for even N, 
( gi)e{x,0 ) = (g i) e (—£,0) and (gi) o (x,0) = (gi) o (—x,0). Expanding Eq. [5(3 however, shows that the limits 
for |x| —> oo and even N are different: (gi) e —> 1, but (gi) 0 —> 1 — 4/N. The intuition behind this result is 
that a strong magnetic held can perfectly align the spins in the even, but not in the odd model. 

Closely related to the nearest-neighbour correlations is the internal energy (i.e. ensemble average of the 
Hamiltonian) per spin 
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-NW lUZ ’ 


(57) 


In general, the calculation yields rather lengthy expressions. However, if H = 0, then U e , 0 = —J(gi)e,o(x = 
0 ,3/). If, on the other hand, J vanishes, then U ej0 = —H(m) et0 (x,y = 0). Using our earlier results of Eq. [47] 
and [55[ 


U e ,o( x = 0,3/) 


( cosh^ 2 —STsinh^ 2 —^ _ , 
l e^(cosh JV / 2 yisinh^/ 2 y) 


U u {x 


0,y) = J 


cosh N 1 y— sinh^ 1 y 

e y (cosh w y+sinhA y) 



if N is even, 
if N is odd, 


(58) 
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H sinh(2a;) cosh^ 2 x ± (—1)^ sinh 


U ey0 (x,y = 0) = - 


N-2 


(59) 


2 cosh^ x ± (—l) w sinh jV a; 

while U u {x , y = 0) = — H tanhx. 

Taking another derivative of In Z with respect to f3 gives us the heat capacity per spin, which measures 
the fluctuations in the energy, 




(60) 


For vanishing H , we can use (dU)/(d(3) = J(dU)/(dy) and Eq.[58l If J = 0, then C = k B /3H 2 x, so the heat 
capacity follows directly from Eq. l50l 


Ce,o(x — 0, y) — 

cosh^/ 2 yisinh^/ 2 y 


COS. 


h N /2 


-2 


y =F sinh^ 2 2 

N-2 • -iN—2 


, N sinh iV/2 ~ 2 (2y) 

^ 2 n / 2 ~ 1 (cosh^/ 2 yitsinh^/ 2 y) 


C M (^ = 0, y) = cosh £ B / +s ^ h v y [ cosh" z y — sinh” ‘y+ 2N .^N y+ ^ y ) 


C'e,o(x-i y — 0) — 

/ 

k B p 2 H 2 


cosh 


■ =F (—1)^ sinl/ 


(-l) Jv iVsinh JV 2 (2a;) 


± 


^ cosh x ± (~1) N sinh * 2 ;v -2 (cosh^ a: ± (-1)* sinh^ x)‘ 

approaching C u {x = 0, y) = k B p 2 H 2 cosh” 2 (a:) in the thermodynamic limit. 


(61) 


if N is even, 
if N is odd, 


(62) 


8. Correlation function 


We can generalize the calculation in the previous section to find the correlation between fc-th nearest 
neighbours. For this purpose we make the spin interactions J in an auxiliary Hamiltonian E dependent 
on the position i. but for simplicity’s sake we drop the magnetic field, 


N 

E(cr) = - ^2 Ji<JiOi+ 1 . 

i=1 


(63) 


Applying the same line of reasoning that led us to Eq. 1351 we can show that the partition function for the 
Hamiltonian E in the case of even N + is 


Ze= J Tr (nQ*) ’ 


(64) 


Ki=l 


where 


Q, = 


/ 0 

e-PJi 

P pJi 


0 


oPJi 


e \ 


(65) 


e -P J i 

eP Ji 0 0 

z~P Ji 0 0 

Vo 0 e _/9Ji eP Ji ) 

plays the role of the transfer matrix of Eq. 1331 The matrices Q; do not commute and therefore we cannot 
simultaneously diagonalize them for computing the trace in Ea. 1641 However, the product QjQj+i commutes 
with Q i+2 Q; +3 . These products are diagonalized as RQjQ i+ iR by the matrix of eigenvectors 


R = 


( 1 1 1 1 

1-1 1-1 
1 1-1-1 
V 1 -1 -1 1 J 


R 
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( 66 ) 
















and the corresponding eigenvalues are 

(RQ i Qi+iR) 11 = 4cosh(/3Ji)cosh(^J i+ i), 

(RQ,Q i+ iR) 22 = 4sinh(^Jj)sinh(/3J i+ i), 

(RQjQ i+ iR) 33 = 4sinh(/?J i )cosh(/3J i+ i), 

(RQ.,Q !+ iR) 44 = 4cosh(/3Ji)sinh(/3J i+ i). 

If N is even, it follows that 

, N N 

Ze, even N = Z N ~ 1 ( COSh(/3Jj) + sillh(/3Jj) + 

^ i=l i=l 

N/2 N/2 . 

cosh(/3J 2 i-i) sinh(/3J 2 j) + sinh(/?J 2i _i) cosh(/3J 2i ) j, 

i=1 i=1 ' 

while for odd N the partition function is half of the unconstrained model’s partition function 

/ n n \ 


(67) 

( 68 ) 

(69) 

(70) 

(71) 


Ze, odd N — 2 Zvl — 2 


N -1 


cosh(/3Ji) + sinh(/3 J,) 


(72) 


For the odd model, we can apply Z a = Z u — Z e . The disconnected correlation function ( g k ) can now be 
computed as 


1 


(»> = v(E 




19 9 


9 




.P k z dJ 1 dJ 2 '" dJ k \ J, = ... =JN=J 


(73) 


with the final result 

rco S h N - fe ysinh fe y+smh iv - fe yco S h fe y±2 1 - JV/2 co S h(2y)sinh iv / 2 - 1 (2y) Ar _ „jj 

-. N/ y 2 - ii iv even, k oaa 

(cosh JV/ ^ y±sinh iV/ ' 2 y) 

coslA^ y sinh fc y+sinh^-* y cosh fc y ± 2 1 ~ JV / 2 s ; 


(3fc)e,o — < 


! sinh iv/2 (2y) 


/ \ _ cosh^ y sinh fc y+sinlA k y cosh fc y 

k \9k/u , 


if IV even, k even, 
if N odd. 


(74) 


(costA / 2 yibsinh^v / 2 y) 
y sinh fc y+sinh^" 18 - 

cosh^ y-|-sinh^ v y 

Taking the limit N —> oo while keeping y and k fixed, all three cases have the same asymptotic value 
liniAr—>oo(fffc}e,o,u = tanl/ y and the correlation length is hence £ = — [ln(tanh y)\. The divergence at y = 0 
can be expressed as a power law in the reduced temperature [2U] 

T r = e~ 2y (75) 

because near T r = 0 


£ w 2 T~ l oc T; 


(76) 


where the critical correlation length exponent satisfies v = 1 in the unconstrained, even and odd model. 


9. Approach to the thermodynamic limit 

It is not surprising that v does not depend on whether we constrain the number of positive spins to even 
or odd values or have no such constraint. We have already pointed out the reason after Eq. [IT] the 
thermodynamic limit at fixed temperature is determined by the leading eigenvalue Ai, and this eigenvalue 
is common to the transfer matrices P and Q. The leading-order correction to the magnetization, however, 
depends on the eigenvalue with the second largest absolute value. If we call this eigenvalue A s , then the 
average magnetization for a chain of length N behaves asymptotically as 

f \ \ N — 1 Q f \ \ 

(to) = si + ( t- ) — ( y— ) + higher order terms, (77) 

V Ai / ox \Xi J 
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Figure 5. The difference between the magnetization (m) and s i for a finite chain of length N, x = 0.5 and 
y — 1. We note that si, defined in Eg. 1431 is the thermodynamic limit of (m). Circles, + and x symbols 
are exact results. The solid and dashed lines are the leading-order approximations of Eq. I78l and l79l The 
constrained models converge much more slowly than their unconstrained counterpart. 


obtained by expanding the logarithm in Eq. 02] and inserting the definition of si from Eq. 03] In the 
unconstrained case we have A s = A2, but for the constrained models one of the other two eigenvalues of the 
matrix Q has a larger absolute value so long as x ^ 0 or y ^ 0. For example, if x and y are both positive, 
then A s = A4 and the leading-order corrections for the unconstrained, even and odd models are 

N 


(m) u - s 1 = -2si 


(78) 


(m)e - s 1 = ~(m)o } N + Si = (ci - Si) 


N 


(79) 


V 

respectively (see Eq. 00] for the definition of Ci). In general, |A41 is considerably larger than |A21- As a 
consequence, the leading-order correction decays much more slowly in the constrained cases than in the 
unconstrained one (Fig. 0. 

The difference in the asymptotic approach to the thermodynamic limit becomes even more apparent if 
we take the limit TV — > 00 while simultaneously J — > 00 (so that T r — > 0) and H — > 0 (so that x —> 0) in such 
a way that the products 


t = NT r , (80) 

h = Nsiuhx (81) 


are constants. In the thermodynamic limit the magnetic field H scales oc TV -1 . We could have alternatively 
defined h = Nx to make this inverse proportionality more apparent, but the definition of Eq.0T]is a little bit 
more convenient when substituting the hyperbolic functions in Eq. 08] and 1391 After applying the formula 
limAr^. 00 (l + z/N) N = e z to Eq. IdOland Eq. 01] we obtain the partition functions for large N, 

Z u = 2 N N/2 t~ N/2 cosh Vh 2 +t 2 , (82) 


N N / 2 t N / 2 (cosh \Jh 2 + t 2 ± coshh) if N is even, 
N N / 2 t~ N / 2 (cosh yjh 2 + t 2 q= sinh ft,) if N is odd. 


(83) 


All thermodynamic quantities can now be derived from the partition function by taking the appropriate 
derivatives, for example (m) = <9(ln Z)/dh. Alternatively, we can also take the thermodynamic limits of 
Eq. 05] 02] and [6T] We tabulate the results in Table [j] 


13 













unconstrained 

even model 

odd model 

(to) 

h tanh y/h 2 -\-t 2 

Vh 2 +t 2 

. h sinh y/h 2 +t 2 - 1-sinh h 

^h2 +t 2 

. h sinh (y/h 2 +t 2 '\ — sinh h 

Vh 2 +t 2 V > -(AT 

cosh V h 2 -\-t 2 +cosh h 

—. h sinh y/h 2 -ft 2 —cosh h 

cosh(Vh 3 +t T )-cosh/t see, 

. h sinh(y/h 2 -\-t 2 )-\-cosh h 

Vh 2 +t 2 ^ n , f AT nAA 

cosh y/h 2 -\-t 2 — sinh h 

cosh(V/i a +t a )+sinh h 

X(h= 0) 

0 

N tanh t 
t 

/tanh(f) r \ 

/coth(|) 1 ^ if y p oven 

[ * + 2 cosh 2 (i) J 

^ tanh t 

( * 2 sinh 2 (|) J 11 iV 1S eVeI1 ’ 

-sshr) if» odd. 

C(h—0) 

k B 0 2 J 2 

n ( tanhi + ^r?r) 

n ( 2tanh (l) + ;^%) 

H (tanhf 

N ( 2 COth ( 2 ) —k%)) lf N 18 eVen 

+ ^t) lf N 18 0dd ' 


Table 1. Thermodynamic properties for TV —> oo and constant t , h (defined in Eg. I80l and l81l) . 



Figure 6. The scaling function of the susceptibility y(/i = 0) in the limit N —»• oo, J3 —> oo with finite 
t = Ne~ 2 P J . The scaling functions for the unconstrained, even and odd models exhibit different behaviour, 
especially if t is small. 


It is instructive to compare these equations with the canonical finite-size scaling forms, for example for 
the susceptibility 

X (h = 0) oc lV/ v f x (N 1 / v t). (84) 

While we find that 7 = v = 1 for all of the cases listed in Table [T] (see also our remark after Eq. I7U1) , the 
scaling functions f x (plotted in Fig. [6]) are fundamentally different. 

10. The probability distribution of the magnetization 

Because of the differences between the unconstrained, even and odd models in Table [T) one may wonder how 
the probability distribution of the magnetization ns cm 121 M differs; after all, (m) and x are essentially 
the mean and variance of this distribution. We repeat here the arguments developed by Antal et al. E5] for 
the unconstrained model with zero magnetic field. We denote the total magnetization by M = ]Th ai and 
the number of domain walls (i.e. boundaries between stretches of contiguous positive and negative spins) by 
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2d ; it must be an even number because of the periodic boundary conditions. The main task is to count the 
number f2(d, M) of configurations with 2d domain walls and magnetization M. Their probability P(d,M) 
in thermal equilibrium with H = 0 will then follow from 

e {N-id)y 

P(d, M ) =--- n(d, M), ( 85 ) 

whose marginal distribution 

P{M) = Y^P{d,M) ( 86 ) 

d 

is the probability distribution we are looking for. 

We can find fl(d,M) with the following combinatorial argument. Let us assume that there are N+ 
positive and 7V_ = N — N + negative spins, and that the first spin is positive. We could for example have 


W 0 + + + W\ - W 2 + + + ... + + W 2 d-1— W2d + ++, 


(87) 


where we marked the positions of the domain walls by w\,..., W 2 d- At the periodic boundary between the 
first and last spin there may not be a domain wall (in the example above there is not), but we will always 
symbolically put wq in front of the chain. We now mentally glue wq, ■ ■ •, W 2 d—i to the next spin in the 
chain (indicated by the braces in Eq. 1871) . In this manner, d negative spins are attached to domain walls, 
whereas the remaining 7V_ — d can be freely placed in the d negative domains. The well-known stars-and-bars 
theorem [26] implies that there are different ways to distribute the negative spins. 

The positive spins require a little more care, because there may not be a positive spin trailing the 
domain wall W 2 d- We can account for this exception by not attaching W 2 d to the following spin. There are 
thus d positive spins attached to wo, u> 2 , ■ ■ ■ , u’2d-2, while the remaining N + — d positive spins can be freely 
distributed into d + 1 segments, namely the positive intervals following wq,W 2 , ■ ■ ■ ,u’ 2 d- According to the 
stars-and-bars theorem, there are ( N + ) different possibilities. 

Because the positive and negative spins are placed independently of each other, the number of 
configurations is simply the product of the binomial coefficients () ■ If wc had started the chain 
with a negative spin, we would have obtained the same expression with the subscripts + and — interchanged, 
so that 


n{d,M) 



( 88 ) 


This expression from Antal et al. [25] is equally valid for the unconstrained, even and odd model. The 
constraints only enter in the permitted values for M whose consequence becomes apparent when we take the 
continuum limit. To this end, we take IV — > oo for a fixed value of d and write m = M/N , so that 

jV 2d 1 


fl(d, to) = 


(1- 


2 \ d ~ 1 

mj 


(89) 


2 2d ~ 2 d\{d-l)\ 

For the time being let us assume that \m\ ^ 1 and thus d ^ 0. We can insert Eq. [55] into Eg. 1551 and 1531 but 
have to bear in mind that changing from the discrete variable M to the continuous variable to generates an 
additional prefactor, which we will call N/Am, 


P(m) 


t N /' 2 A M Z(l — m 2 ) 


oo 


E 


t 2d ( 1 — m 2 ) d 
2 2d d\(d — 1)!' 


(90) 


Here Am is the step size between consecutive values of M (i.e. Am = 2 in the unconstrained, Am = 4 in 
the even and odd model) and t is defined in Eq. [Ml As noticed in Ref. [25] , the infinite series in Eq. [90] can 
be expressed in terms of a modified Bessel function of the first kind thanks to the identity [27] 

00 -2d+l 

A(*) = E^XT3F7TT-Tu (9!) 


<f=0 


2 2d+l d\(d + 1)! 
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and therefore 


P(m) = 


2N N / 2 Ii (ty/T — ra 2 ) 

f ZVT~ 


(92) 


t N / 2 ~ 1 A M Zy/L — m* 

for \m\ < 1. 

At the boundaries of this interval (i.e. |m| = 1) there are contributions proportional to Dirac delta 
functions. These singularities arise because \m\ = 1 implies d = 0, leaving the denominator in Eq. [89] 
undetermined. For the unconstrained model as well as the even and odd model with even N, the 
proportionality constants in front of the delta functions can be computed based on the observation that 
P{m) must be normalized and symmetric about m = 0. For the even model with odd N, a discrete 
magnetization M = —N is permitted, but M = N is not, so that a delta function can only appear at 
m = — 1, but not m = 1. Conversely, the odd model with odd N can only have a singular contribution at 
m = 1, but not at m = — 1. 

The probability contained in the regular part of the distribution given by Eq. [22] follows from the integral 

r 1 h 


J ) 


i-i \/l — m 2 t 

and, upon inserting the partition functions of Eg . 1821 and 1831 with h = 0, we obtain 

th ( t\/l - m 2 ) S(m - 1) + <5(m + 1) 


dm = 


4sinh (f/2) 


P u {m) = 


P e (m) = 


2y/l — m 2 cosh t 

tli (ty/l—rri 2 ^ 

4\/l —rn' 2 cosh 2 ( t /2) 
tli(ty/l-m 2 ) ^ 2 8 ( 
2\/l —m' 2 cosht 


cosht 

. 8(m— l)+<5(m+l) 
cosh 2 (i/2) 


if N is even, 


t+i) 


(93) 


(94) 


(95) 


cosh t 


if N is odd, 


tli (ty/l—m' 2 ') 


if N is even, 


Po(m) = { 4Vl- ^ sinh 2 (t/2) “ ’ (96) 

2S(m-l) iiNisodd 
2CT^cosh t + cosh t 11 iv lb OUU - 

One noteworthy detail is that the delta functions peak exactly at the boundaries of the interval [—1,1]. 
So long as the integral of the delta function over the entire real line equals 1, it is a matter of definition how 
much weight is assigned to the left and right of the interval boundaries. We have adopted here the symmetric 
convention J 0 °° 5(x)dx = 1/2 which applies, for instance, if the delta function is the limit of narrowing zero- 
centred Gaussians. Other conventions are possible; for example Ref. [2^ implicitly uses / Q °° 8(x)dx = 1 
which changes the prefactors in front of the delta functions in Eq. 19411961 With our definition of the delta 
function and the integral 

i 2 ) , 2 (sinh t 


rl m 2 /i (ty/1 — ' 


-dm 


■ y/l — rn 2 t \ t 

we can indeed retrieve the susceptibility x in Table [I] 


- 1 


(97) 


11. Discussion 


Combining the results above, we can now analytically solve the stochastic synchronous and asynchronous 
voter models introduced in Sec. [2] and 2] With Eq. [5] we can translate m and gi of the Ising model into 
correlations between the opinions of nearest and next-nearest neighbours, 


<">> = ¥ 




Wit+z+l 


(98) 


(5i) = 


N 


y/ wicuj+2 


(99) 
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Figure 7. (a) Mean nearest-neighbour correlation (m) and second-nearest neighbour correlation (gi) in 
the stochastic voter model. Black curves and symbols are for synchronous, red for asynchronous updates. 
Analytic predictions (Eq. Il02l 110411106111081) are shown as solid and dashed curves. The results of Monte 
Carlo simulations for a chain of length N = 100 are shown as circles and squares, (b) The same for the 
variances of m and g\ fEa. 1103111051110711109|) . 


where the second equation follows from Eq. l5ll and u;? +1 = 1. The variances of m and g\ are proportional to 
the second partial derivatives of In Z with respect to either x or y, thus 

IVvar(m) = 4 (100) 

P 

Nvai( gi )= (101) 

where the derivative in the last equation has to be evaluated at x = 0 for synchronous and y = 0 for 
asynchronous updates. For synchronous updates, we have in fact evaluated this derivative already in Eq. RTT1 
because in this case N vai(gi) = C /(fcs/3 2 J 2 ). The corresponding calculation for asynchronous updates can 
be performed by differentiating the partition functions in Eq. [40] and 01] 

Inserting Eq. [9] into Eq. 0S] 0H1 EH and [61] we obtain for the synchronous voter model in the 
thermodynamic limit 


1 1 if AT is odd and either p + = 0 or = 1, 

lim (m) = \ 

N—too ' 1 

1 0 otherwise, 

(102) 

toJATvarfm)] = 2 ^_, 

(103) 

i™ 

N-i- OO' 1 + 2y/p + P- 

(104) 

lim [Vvarl'oi'l] = -- VP+P - 

oo L (1 + 2 v ^pT) 2 ’ 

(105) 
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where, as before, p- = 1 — p+. For asynchronous updates the corresponding results are 

y/P+~ y/pZ 


lim (m) = 

N ^oo y/P+- 


y/P—’ 

4 y/P+P- 


lim \N var(m)l = „ _, 

N-¥0 O 1 + 2y/p+p_ 


N 


lim (#1) = 


N- 


lim [A r var(g 1 )] = 


1 ~ ^y/P+PZ 
1 + 2 y/p+p- ’ 

^y/P+P-iX ~ y/P+P -) 


( 106 ) 


(107) 


(108) 


(109) 


1 + 4 v ^pT(l + y/p^pZ) ' 

We plot Eg. 1102111091 in Fig. [3 As a numerical confirmation we include the results of Monte Carlo simulations 
in the same graphs. The numerical and analytic results are in excellent agreement. 

Comparing the synchronous with the asynchronous case, we notice that the thermodynamic limits of the 
nearest-neighbour correlations (m) differ significantly. While for synchronous updates nearest neighbours 
are typically uncorrelated, asynchronous updates build up non-zero correlations. Interestingly, the mean 
second-nearest neighbour correlations ( g\) are identical for both update rules. However, the variances differ 
between the rules: var(m) is larger for synchronous updates, whereas var(gi) is larger for asynchronous 
updates. 

It is in principle possible to extend the calculations to correlations between more distant neighbours too. 
For example, (X) uiiUJi+X) can be obtained by formally introducing a three-spin interaction strength K in the 
Hamiltonian so that E(<t) = —K X)j <r i <Tj + iCi+2 — J X)i l — H X) a i- The mean third-nearest neighbour 

correlation follows from differentiating the partition function Z with respect to K and subsequently setting 
K = 0 as well as H = 0 for synchronous, J = 0 for asynchronous updates. Unfortunately, the transfer matrix 
method developed in Sec. [5] does not easily generalize to arbitrary fc-spin interactions [28], but calculating 
correlations in the voter model from Ising-like Hamiltonians is an intriguing possibility for future research. 


12. Conclusion 

We have studied two variants of the one-dimensional Ising model: in the first variant the number of positive 
spins is constrained to an even number; in the second model this number must be odd. We have motivated 
both models by mapping them to a model of opinion dynamics with either synchronous or asynchronous 
updates. If the temperature and magnetic field are held constant, the thermodynamic limits of the even 
and odd Ising models are the same as the limit of the unconstrained model. However, by simultaneously 
increasing the chain length and lowering the temperature and magnetic field, we have shown that the scaling 
functions for the even, odd and unconstrained models differ. The mapping from the Ising model has allowed 
us to obtain explicit formulae for correlations between nearest and next-nearest neighbours in the voter 
model. 

We can generalize the problem posed in this paper to higher dimensions or complex networks by 
associating spins with the links and enforce an even number of positive spins on every cycle in the graph. In 
other words, only balanced signed graphs [29j are permitted. Assigning opinions to the nodes and mapping 
them to spins on the links as in Eq. [5] will naturally generate such graphs from the voter model. It is a 
fascinating question how this changes the thermodynamic limits compared to the unconstrained model. 
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